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Abstract. We study the stability of charged rotating black holes in a consistent 
truncation of Type II B Supergravity on AdS^ x that degenerate to extremal black 
holes with zero entropy. These black holes have scaling properties between charge and 
angular momentum similar to those of Fermi surface-like operators in a subsector of 
Af = 4: SYM. By solving the equation of motion for a massless scalar field in this 
background, using matched asymptotic expansion followed by a numerical solution 
scheme, we are able to compute its Quasi-Normal modes, and analyze it's regime 
of (in) stability. We find that the black hole is unstable when its angular velocity 
with respect to the horizon exceeds 1 (in units of 1/lAds)- A study of the relevant 
thermodynamic Hessian reveals a local thermodynamic instability which occurs at the 
same region of parameter space. We comment on the endpoints of this instability. 
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1. Introduction 

The AdS/CFT correspondence, first formulated by Maldacena [1], lias been used to 
study the microscopic degrees of freedom of black holes, which are known to possess 
a large entropy. Indeed, significant progress was made in understanding detailed 
properties of supersymmetric (SUSY) black holes, via their CFT duals. However, the 
case of little or no supersymmetry is much less understood, and it is not always possible 
to pinpoint the microscopic duals of black holes. 

The most studied example of AdS/CFT correspondence is the duality between 
type II B string theory on AdS^ x S^ and Af = A SU (N) SYM, and it is also the one 
we have in mind in this paper. We will be interested in black holes that have angular 
momentum in a two-plane, which we denote by J^, and two equal charges Qi = Q2 = Q- 
Both the angular momentum in the transverse two-plane and the third U{1) charge are 
set to zero. The extremal limit of these black holes has zero entropy. For a specific one 
dimensonal set of Q and the extremeal limit has 1/16 of SUSY preserved, but other 
than this measure zero set, no SUSY is restored. 

The fact that these black holes have zero entropy, and the details of the scaling 
relations between angular momentm and charge (in the large charges limit) makes them 
ideal candidates for Fermi surfaces in the dual A/" = 4 superconformal field theory. 
Indeed, in [2] we made a conjecture regarding a specific CFT dual configuration, relying 
on the structure of closed subsectors - the PSU{1, 1|2) sector - of the dual theory. In 
order to support the conjecture, we studied the spectrum of scalar perturbations, and it 
was shown to match our expectations when the black holes are fast rotating and close 
to extremality. In some more detail, the appropriate dual was conjectured to be an 
effectively 1+1 chiral Fermi surface, and the spectrum of the quasi- normal modes of the 
scalar above matched the expected dimension from the CFT. 

Further understanding of the proposed duality requires a better understanding of 
the phase diagram on both sides. This paper presents an analysis of the stability, 
or more precisely instability, of these black hole backgrounds, and shows that they 
are prone to two types of instabilities - superradiant instability and thermodynamic 
instability - in part of the parameter space of and Q. We work in the double limit of 
large charge/angular momentum and vanishing temperature, taking first the charge to 
infinity and only then the temperature to zero. In addition, we discuss the near-horizon 
limit of these black holes. The extremal black hole has a pinching AdS^ orbifold near 
horizon, as discussed in [3l H] , where it was shown to appear in the extreme Kerr black 
hole in AdS^, when one of the angular momenta is set to zero. 

To check for the existence of a superradiant instability we add a massless scalar to 
the black hole background and perform a linear stability analysis. This scalar, which 
is the lOD axion-dilaton, can be added within a consistent truncation of the type IIB 
theory on S^. We compute the Quasi-Normal Modes (QNMs) of the scalar in the black 
hole background numerically for different values of parameters in the above-mentioned 
limit. 
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To the best of our knowledge, linear stability analysis has not been performed 
for these black holes. However, there have been several studies aimed at computing the 
scalar field QNMs of rotating black holes in AdS. Analytical results have been obtained, 
but only for slow- rotating small black holes (as compared to the AdS radius), where one 
can neglect the effect of the black hole far away from it. The analysis was carried out in 
-D = 4 [5] and in D = 5 |6]. The fast rotating large black hole regime has proven more 
difficult to explore, even numerically [7]. 

At the linear level, the equation of motion for the added scalar is that of a minimally 
coupled massless scalar field. The equation is separable, reducing it from a single PDE 
to a pair of coupled ODEs, a radial and an angular one. For a generic choice of the black 
hole's parameters, these are both Heun equations, i.e., second order regular differential 
equations having four regular singular points. In the fast-rotation limit, two of the 
singularities in the angular equation merge. Then, using a matched asymptotic analysis 
developed in [HI |2], with some modifications, we are able to solve this equation. The 
solution yields the separation constant in terms of the complex frequency of the scalar. 
Finally, we use the Continued Fraction Method (CFM) to solve the radial equation 
numerically and find eigenfrequencies. 

It was shown in [9] that for Kerr- AdS black holes there can be no superradiance in 
AdS^ if the angular velocity at the horizon, denoted Q^, is smaller than one. For > 1, 
on the other hand, one expects the appearance of supperradiance, which in an AdS 
background renders the black hole unstable. Our expectation is that this remains true in 
the large charge/angular momentum limit, so that the black hole becomes superradiant 
at > 1. However, in the large charge/large angular momentum limit of our black 
holes, the angular momentum is parametrically larger than the charge and i7 — )■ 1 for 
all black holes in this limit. We therefore need to examine more closely from which side 
it approaches 1, ClS Sb function of the subleading charge. 

Parameterizing by a the ratio between the angular momentum and the charge 
squared, a = we have, close enough to extremality, < \ when a < 2 and 

> 1 when a > 2. For a = 2, fi^ = 1 and the extremal black hole satisfies a SUSY 
BPS bound, i.e., the SUSY BPS black holes divides the parameter space into two regions 
with very distinct features, with expected stability on one side and instability on the 
other. Indeed, our numerical analysis verifies this picture and we find no superradiant 
instability for a < 2. For a > 2 instability is possible only for modes with very 
large values of angular momentum, ~ ^J~J^J~N^ of the scalar. For such values, the 
perturbative expansion we used to obtain the solution is under lesser control. However, 
the application of the solution in this case displays superradiance accompanied by an 
instability. 

Next we turn to the issue of a thermodynamic instability. To determine wether 
the black hole is thermodynamically stable we have analyzed the eigenvalues of the 
thermodynamic Hessian. For fi^ < 1 we find that the black hole is thermodynamically 
stable, i.e. all eigenvalues of the Hessian are positive . For fi^ > 1 there is one negative 
eigenvalue corresponding to the unstable direction, characterized by an equal change in 
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energy and angular momentum and no change in the charge. We thus find that both 
thermodynamic and superradiant instabihties are present for black holes with fi^ > 1, 
and no instability is observed for black holes with < 1. This is reminiscent of the 
Gubser-Mitra conjecture, but for rotating black holes [10]. We will return to this point 
in the summary and discussion section. 

The outline of this paper is as follows. In section [2] we briefiy review the stability 
problem for black holes, specifically focusing on the phenomenon of superradiance and 
superradiant instability. Section |3] is dedicated to describing the black hole background 
and examining its extremal limit, followed by a discussion of the near-horizon limit 
of such geometries and a short description of the CFT dual. In section H] we solve 
the angular equation using an asymptotic matching technique, and set up a numerical 
procedure for solving the radial equation using the Continued Fraction Method. The 
results of the computation are presented in section [5l where we compute several 
Quasinormal modes using the method described here. In section [6] we analyze the 
thermodynamic instabilities of the rapidly rotating black hole. Finally, the results, 
along with possible end-points for the instabilities, are discussed in section [3 

2. Black hole stability 

The stability of black holes is the subject of extensive research. Usually, the question of 
stability is divided into two parts: classical stability and quantum stability. The former 
can be studied in general relativity, usually in linear perturbation theory with higher 
orders for zero modes. The latter can be understood using semiclassical methods such 
as black hole thermodynamics. In this paper, we discuss both types of (in) stability, with 
emphasis on classical stability. 

2.1. Classical Stability 

One usually studies the classical stability of black hole backgrounds by using the 
linearized equations of motion, within the low-energy effective action, for whatever fields 
are present in the theory. While some works have attempted to solve the full non-linear 
stability problem and find the final stable configuration, this is usually a very difficult 
task which we shall not discuss here. 

In order to analyze the linear stability of black holes, one first writes the linearized 
equations of motion around the black hole background for the different fields present in 
the (super)gravity background. Typically, this system is not tractable analytically, or 
numerically (with any reasonable amount of effort). However, in some cases, the various 
equations decouple, so that the stability can be assessed for each field (or combination 
of fields) separately. As black holes are time independent, one can then look for the 
eigenmodes of time translation for each field on this background. Since we are dealing 
with a dissipative system, as all waves have to be incoming at the event horizon, these 
are not normal modes, and are rather called Quasi- Normal Modes (QNMs for short). 
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having complex eigenvalues under the time translation symmetry (for a review of QNMs 
and their role in black hole physics see [H]). Stability follows directly from the sign of 
the imaginary part of these eigenvalues - in our convention, when the sign is positive, 
the perturbation decays, and when it is negative, the perturbation grows, signaling an 
instability. 

One of the most well-studied instabilities of this type is the superradiant instability, 
occurring in charged or rotating black holes. The black holes discussed in this paper are 
both charged and rotating, and hence this phenomenon merits special attention. 

2.1.1. Superradiant Instability The phenomenon of superradiance allows the extraction 
of energy from a black hole by scattering specific modes, whose reflection coefficient 
from the event horizon is larger than one. It was first discovered by Zel'dovich [T2] . 
who noticed that a conducting cylinder rotating about its axis with angular velocity 
Vt can amplify electromagnetic modes impinging on it which satisfy uj — mQ < 0, m 
being the angular quantum number and u the eigenmode's frequency. This process thus 
extracts rotational energy from the cylinder. There have been many results identifying 
superradiance in different types of black holes: 

• Bekenstein [13] has shown that if the area theorem holds then modes satisfying 
Rew < m^fl^ will lead to the extraction of energy from the black hole. Here, 

is the conserved quantum number of the azimuthal angle 0, with respect to a 
non-rotating frame at infinity. Also, we denote by the angular velocity of the 
black hole at the event horizon, defined as the coefficient of in the killing vector 
which becomes null on the horizon, I = + il^-^. It was then shown by direct 
computation that indeed scalar field modes satisfying Rew < m^Q^ impinging on 
the black hole are superradiant [H]. The same can be shown for gravitational and 
electromagnetic perturbations [T5] . 

• Consequentially, it can be shown that given boundary conditions that correspond 
to refiecting the superradiant modes back onto the black hole, say by placing it in 
a box, it will grow unstable. This is referred to in the literature as the 'Black Hole 
Bomb' [16j. One thus expects black holes in an AdS background with large enough 
angular velocity to be unstable due to superradiance. 

• Indeed, one can also show by explicit computation that in small slowly-rotating 
Kerr-AdS black holes in Ads4^ a scalar field mode is unstable iff it is superradiant 
[5]. This result has been obtained for gravitational perturbations as well |17j . 
These solutions have been obtained by solving the angular equation perturbatively 
about the flat-space solutions, using the slow rotation and small size of the black 
holes. Similar analysis has been performed for charged rotating AdS^ black holes 
[6]. When these conditions are not satisfled, an analytic solution is typically not 
available, yet several numerical methods have been used to flnd QNMs, displaying 
instabilities in some cases (e.g. [7]). 
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In this paper we explore black holes in an asymptotically AdS^ spacetime, which 
are rotating with angular momentum scaled to oo (around one of the rotation planes). 
The size of these black holes is comparable to lads- For such black holes, there are 
no known analytical solutions for the QNMs. There has been an attempt to solve the 
linearized equation of motion numerically for the Kerr- Ads black hole, and even in this 
case one encounters problems in the fast-rotation regime Here, we apply a technique 
introduced in [8j and refined in |2j to find semi-analytically the QNMs for fast-rotating 
black holes in AdS^. 

2.2. Semiclassical Instabilities 

Black hole thermodynamics goes back to Bekenstein's and Hawking's seminal works 
[T8| [To] . The thermodynamic effects in black hole backgrounds can be studied using 
semiclassical calculations. There are two types of instabilities which are discussed in 
this context: 

• A gravitational solution may not be a local minimum of the free energy, in 
which case it is thermodynamically unstable. The prime example of this is 
the Schwarzschild black hole. A generalization of this, for multi-charged black 
holes in various dimensions has also been studied (e.g. [20]). Specifically, such 
thermodynamic instabilities manifest themselves as non-positive modes of the 
thermodynamic Hessian 

H - m 

where the extensive variables describing the black holes, for example 

the total energy, Q, the electric charge and J, the angular momentum. While this 
computation is straightforward in simple backgrounds, such as Schwarzschild and 
Kerr, it becomes more complex as charges and angular momenta are added in more 
complicated gravity theories. 

• Global instabilities: these occur when at least two gravitational backgrounds with 
the same charges exist, and are local minima of the free energy. In this case, the 
relevant free energy may be lower for one of them compared to the other, rendering 
it the stable configuration, with the other being metastable. If, in some region of 
the parameter space, these roles are reversed, a phase transition will occur. The 
most well known example of this phenomenon is the Hawking-Page phase transition 
|21j which has been given a beautiful interpenetration within the framework of 
AdS/CFT, relating it to the confining/deconfining phase transition in the dual 
gauge theory [2^ . 
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3. The Black Hole 

3.1. Description of the black hole solution 

The notations in this paper follow closely those of [23]. We work within the consistent 
truncation of type II B Supergravity on AdS^ x described in [21]. The field content 
comprises of the metric, two neutral scalars and three abelian U{1) vector fields. It 
can also be consistently coupled to a dilaton and an axion, which are massless, as was 
discussed in [21] and later used in p]. The bosonic part of the supergravity Lagrangian 
is given by 

2 3 



1 



a=l i=l ^ ^ 



+ 



24 



\ uvpu -pi -pj Ak /r,\ 



Here, lads = / = - is the AdS radius, A^ are the three U{1) gauge fields, and Xj 
are the three uncharged scalars, constrained by X1X2X3 = 1. These scalars may be 
parameterized by 

Xi = e-75^^-75<^^ X2 = e"^*^+^'^^ X3 = e75'^^ (3) 

We discuss the black hole solutions that first appeared in |25J, while following the 
notations of [23] which contains a generalization of these solutions. The solutions are 
parameterized by 5i, ^3, m, a, 6, which map to Qi = Q2, Q3, E, and Jr. Motivated 
by dual field-theoretic considerations first presented in [2], which we describe briefly in 
section |3^ we are interested in solutions for which Q3 = 0, = and Qi = Q2. This 
amounts to setting 6 = and ^3 = 0. The solution then reduces to the following metric 

,,2 _ rrV'sf (2 2^ [^,^\_ Xjdt-y^daf 

Y \dt + [x^ + 2ms'i) daf' 99,9! , 



(a;2 + y2) Hi J 
were we deflne si = sinh{6i) (and ci = cosh {61) for later use). The functions used in 
the metric are 

X = -2m + (a^ + x^) + g"^ (a^ + 2msl + x^) {2msl + x^) (5) 
Y ={a'- y') (1 - gV) (6) 
2msl 
+ y'-^ 

The linear transformation to the coordinate frame that is non rotating at inflnity 
is given by 

t =^^- (8) 

X = -i^, (10) 

a 
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with = 1 — g^a"^- Finally, it is occasionally convenient to use an additional 
transformation in the non-rotating frame 

2 2^2 /I 1 N 

X = r ms^ [11) 

3 

y2 = (12) 

The event horizon is now located at r = ro, the largest root of 

A(r) = (13) 

The angular velocity around the direction at the event horizon, with respect to 
the asymptotically non-rotating AdS space, is given by 

Conversion to natural field-theoretic units, where Qi = Q2 = Q and Q3 are the cartans 
of 5*0(6) and J^, are half- integer quantized, is performed using the relation [26] 

if = T- (1^' 

In terms of these units, the black hole's charges are [23] 

vr 2ma(l + g^) 2^3 ma(l + g^) 

g.-Q2^g-4^-iVV^ (17) 
vr m[{2{g^a^ + + 1) + g'a\^a - 2))sl + + 2] 

^ = 4g: ^^^^ 



^2 2^rg - |msf(a2 + |msi 



N'g'n^ (19) 

ro (drA) \Jrl- fmsi^ 
A-K {a? + rg + Imsi^) (rg — Imsi^) 

Here, ro is again the radius at which the event horizon is located. In the rest of this 
paper, we set / = 1/(7 = 1. One can restore these factors using dimensional analysis. 

3.2. The Near Extreme and the Fast Rotation Limits 

First, we briefly analyze the extremal limit of the black hole. In order to find the location 
of the horizon we solve A(r) = in f[T3|) . There are two solutions for A = 0, given by 
X = or = 0: 

x- = ^:rl = ^^ (21) 
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X = 0: rl = ^ i ^ (22) 



rl = ^ ^ ^ ^ ^ < (23) 

There are two horizons, located at ri and r2 and the extremal limit is obtained when 
ri = r2 in which case we also find that the entropy is zero. Using these definitions we 
write 

X = (r^-rl){r^-rl) = (x' + {rl-rl)){x^ + {rl-rl)) = (z - z^) (z - Zi) , (24) 

where in the last line we denoted 2; = x^, so that Z2 = r2 — rf and zi = r| — r^. When 
r2 < ri we get a naked singularity. Thus, for the non-singular black holes we can use 
Z2 as a measure of the distance of the black hole from an extremal one. It will be 
convenient to use a parametrization in terms of a = j^^^^, a and Z2- In terms of these 
parameters, we have 

{a^ + Z2) {I + Z2) , 2a2 a3 + 2az2 
^ = 2 + ^ + ^^ ^25) 

(26) 



a 
a 



Aa + {a^ + 1 + Z2) a 

Zl — . [Zl) 

a 

In this parametrization, the near extremal limit for this black hole is thus Z2 — 0. 
We have studied this limit in [2], where we have found some analytic solutions for the 
QNMs, which are similar to those of a 1 + 1 CFT. In the class of large rotating black holes 
that we consider here, there is an additional limit which facilitates solving the equation 
of motion, namely the fast rotation limit. For black holes with angular momentum 
scaled to infinity, the distance of the rotation parameter a from its maximal value 1 can 
be used to obtain some analytical results. In this paper we focus on the large rotation 
limit, and thus Z2 ^ (1 — a). 

The BPS limit of the black hole is 

El = J^ + J^ + Qi + Q2 + Q3 = + 2Qi (28) 

and it is satisfied at 2:2 = and a = 2. We expect superradiance, and since the spacetime 
is asymptotically AdS also an instability, if fi^ > 1. Rewriting fi^ in terms of a and Z2 
we get 

= + " + ^"^^ = 1 + - " + ""^^ - + O ((g - 1)^) . (29) 
a(2 + aa) + az2 2 + a + az2 

This expression shows that to leading order in (1 — a) we have = 1 independently 
of both a and Z2- However, the next-to-leading order correction shows that for fixed 
a there is a threshold value for the off-extremality parameter Z2 = 1 — 2/a, such that 
Q(f, < 1 for Z2 > Z2 and i7<^ > 1 for Z2 < Z2. Since Z2 is always positive, we can see 
that for a < 2 we have fi^ < 1 for all values of Z2. As the appearance of superradiant 
instability is expected only for fi^ > 1, in the limit a — t- 1 it should appear only for 
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a > 2 and sufficiently low Z2. Note that in the extremal limit, i.e when we set Z2 to zero 
ffist, the crossing point from > 1 to fi^ < 1 is precisely at a = 2 which corresponds 
to the BPS black hole. 

It is useful to present the temperature and entropy in terms of Z2 and a: 

^27r^(a'^ + ^ + Z2) , , 

S = ^ ^ " (30) 

y _ y/z^ (4a + g + g^g + 2^2:2) 
~ Aan + 27ia (a^ + Z2) ' 
Let us express the appearance of superradiance in terms of the temperature. To 
leading order in 1 — a we have 

^ + 0[l-a]. (32) 
vr 

There is then a critical temperature Tc = ^ + 0(1 — a) above which is smaller than 
one for any a and we expect no superradiant instability to develop, independently of 
the choice of a. 



3.3. The Near-Horizon Limit 

It is usually useful to study the near-horizon limit of black holes in order to understand 
their properties. Such studies, mainly of extremal or near-extremal black holes, have 
resulted in better understanding of the black hole's microscopic degrees of freedom. The 
issue of instability can also be studied using the near horizon geometry, which is often 
simpler than the full geometry [27j. Thus, it would be beneficial to find the near horizon 
geometry of the black holes considered here. There is, however, an obstacle to defining 
this limit globally for the near extremal black hole, as we discuss below. 

Examining ( 130|) we see that the entropy, and thus the horizon area of the black hole, 
approaches zero together with the off-extremality parameter Z2. The extremal metric is 
therefore singular. In addition, computing the Ricci scalar, one can see that at x = 
there is a ring of singularities at y = 0, where the ip circle pinches to a point. This ring 
is hidden behind the horizon for the non-extremal black holes, and touches the horizon 
when Zo = 0. Below we present the extremal near- horizon limit for this family of black 
holes, lil 



The Extremal Near-Horizon Geometry a; — ?■ 0, fixed y: To take this limit carefully we 
ffist define x = ep, taking e — )■ in the end. Applying this to the functions appearing 
in the metric we get: 

X ^ (^l + a2 + 4^)eV (33) 
H,-^l + —^. (34) 
I The analysis of this section was carried together with Joan Simon, and wiU be discussed further in 

EHl 
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If we naively take e — )■ we will set the time coordinate of our solution to zero. Instead, 
we take t = K We will also define (p = t + and x = ^X- The metric then takes the 
form: 

ds' = H!h(, f , , +p'df-(l + a^ + A^) ?^p'dA+y^+^].m 

This metric has a pinching AdS^ orbifold factor, as defined in [29], which is an AdS^ 
except that the periodicity of the spatial circle x scales to zero in the this limit. Such 
geometries have been studied in the context of Extremal Vanishing Horizon (EVH) black 
holes in AdS^ [H [3] , and recently in [30] . It would be interesting to see whether one can 
use this generic structure, along with the proposed dual, to learn more of the extremal 
black hole's properties and test the EVH/CFT duality in detail. 
Subsequent y — 0: Taking now the ?/ — )■ limit yields: 

We see that the radius of curvature of the AdS^,, and actually of the entire space, 
approaches zero when ?/ — > 0. Hence ?/ = is a naked singularity in the extremal case. 

The loci ?/ = is also singular in the following sense: the distance to x = at y 7^ 
is infinite in the extremal limit, whereas the distance to x = at y = is finite. The 
situation is exasperated by the fact that the distance to ?/ = along the x = horizon 
is finite even at extremality. 

The computation is straightforward. The metric component in the radial direction 
is given by 

Now, taking the limit x — )■ with finite y, we obtain 

, 2a y' 1 

We see that the distance to the horizon, located at x = 0, along the radial direction is 
infinite. This is the familiar 'throat' of extremal rotating black holes. 

On the other hand, approaching the x = limit at fixed ?/ = 0, the metric 
component in the radial direction is now given by 

2a V'^ 1 1 



axy x2 + 1 + Amsl + a^ ^ x^/s ' ^^^^ 
Performing the integration from finite x to x = now yields a finite proper distance. 



3.4- Field Theory Dual 

For completeness, we would like to present the motivation for studying this class of 
black holes, which is rooted in the AdS/CFT duality In [2] a conjecture was made 
regarding the field theory dual of these black holes. The dual conformal field theory 
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for this supergravity theory is J\f = 4 SU{N) Super Yang-Mills (SYM). In this theory, 
there are several known sectors which are closed under the dilatation operator, which 
is the radial quantization equivalent of the hamiltonian [31]. Thus, a state composed of 
the operators in such a sector will only mix with other states within the sector. One of 
these sectors, named the PS'?7(1,1|2) sector, contains only partons with the following 
relations between the charges 

Ao = 2Jl + Qi + Q2 + Q3 = 2Jr + Q1 + Q2-Q3 (40) 
where Aq is the classical scaling dimension, and Jr are the SU (2) x SU (2) quantum 
numbers, and Qi, Q2 and Q3 are the SU{4:) i?-charges spanning the U{lY Cartan 
subalgebra of S0{6) = S'[/(4), which are, of course, the U{1)^ charges discussed before 
in the AdS^ x supergravity dual. 

The Af = 4 SU{N) SYM theory contains a single supermultiplet in the adjoint of 
the SU{N) gauge symmetry. This multiplet includes a gauge field A^, a fermion ipa 
and its conjugate ip"', in the fundamental and anti- fundamental representations of the 
S'f/(4) R-symmetry, respectively, and a scalar cp^ab] in the 6 of SU{4:). The lower indices 
on these fields are fundamental indices of the SU{4) R symmetry, while upper indices 
are anti-fundamental ones. 

The constraints fHOj) retain only four types of partons from the full set 

= D^^^t/'uV'^") = V1-4)\. (41) 

Here, D^^ is the covariant derivative D^a (exchanging the 5*0(4) of Euclidean rotation 
with SU{2) X SU{2) notation) with a = a = 1, and n indicates the number of derivatives 
operating on the fields. 

Each of these fields also sits in the adjoint of the gauged SU{N) so it has an 
additional adjoint index which we have so far suppressed. We will sometimes wish to 
write this index explicitly, e.g. 0"^.^ where j = l...dimG = N"^ — 1. One can build gauge 
invariant quantities in the following way. Given an operator \& in the adjoint of SU{N), 
define Jdet[^] as 

9 

Jdetm = e'^i'^^-'^«^(,,)^(a.)...^(a,)* = Yl ^W^" (9 = dimG). (42) 

a=l 

The charges of the fermions in this sector are given in the following table. 

Table 1. Quantum numbers for the fermionic content of the PSU(1,1\2) subsector, jp 
refers to the angular momentum of the SU{2)p automorphism. 





jL 




Qi 


Q2 


Q3 


jp 




(n + l)/2 


n/2 


1/2 


1/2 


1/2 


1/2 




n/2 


(n + l)/2 


1/2 


1/2 


-1/2 


-1/2 



The pair of fermions in UT] enjoy a further SU{2) automorphism, under which ip^"^^ 
and '?/'*■"'' transform as a doublet, and we will denote them \E'^"''* and '^'^\ respectively, to 
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make the SU{2) action manifest. This symmetry can be employed in order to build 
states which cannot mix with any other state under the dilatation operator. One 
can build such states which have the same quantum numbers as the black holes we 
examine. This has been done in detail in [2], where states resembling a Fermi surface 
have been constructed. These satisfy the field theory BPS bound at zero coupling, and 
are conjectured to be duals of the extremal BPS black holes within this family, i.e. those 
having a = 2. 

Using the notations described before, the basic operator, corresponding to the a = 2 
extremal black hole, is defined as 



Here, SymW stands for a symmetrization of the operator with respect to the SU{2) 
doublet indices {1,2}, placing the operator in the highest SU{2) pseudo-spin state, 
with = 2KN'^{2KN'^ + 1) and Jz = 0. Higher values of alpha can be obtained by 
removing some of the operators in 0'^^\ thus puncturing holes in the Fermi surface. 
While there may be several ways to obtain the same set of charges, the entropy is still 
0(1) and not 0(A^^), so it cannot be seen on the gravity side. 

4. Scalar Perturbations in the Black Hole Background 

In this paper we are interested in near extremal, rapidly rotating black holes such that 
Sa = (1 — a^) ^ y/z2 which results in a black hole with temperature close to zero, 
but non-zero entropy. We wish to examine the possibility of superradiant instability 
in these black holes and extract information about the QNM of perturbations on this 
background. 

In order to be able to interpret the resulting perturbation spectrum in the dual 
CFT, we would like for the perturbing field to be part of the consistent truncation of 
the full type II B theory. While there are two scalars in the Lagrangian (|2]), they are 
both massive, and their equation of motion is not separable. However, as shown in 
[32], it is possible to use the symmetries of the ten-dimensional theory to extend the 
consistent truncation considered thus far, allowing us to include the axion and dilaton, 
which transform under the SL{2, M) symmetry of II B. 

At the linearized level, the equation of motion for the dilaton decouples and turns 
out to be that of a minimally coupled massless uncharged scalar. Thus, this is the 
simplest field we can consistently turn on to check the existence of an instability. The 
equation of motion is given by 



VK-l 2K-1 




(43) 



n=0 m=K 




(44) 



Plugging in the ansatz 



^{x,y,T,(f),ip) = e 



"»^+^'r-^^S{y)R{x), 



(45) 
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the equation of motion can be separated into two ODEs. This separation of variables 
introduces a new parameter, A, in addition to m^, and u. We will carry out the 
analysis of stability in the limit limit — with ^ ^/z2■ Some aspects of the 
analysis are similar to that of j2j. However, there the opposite limit was considered, 
which in some sense 'reverses' the roles of the angular and radial equations. 

The outline of this section is the following. In section 14.11 we examine the angular 
equation, i.e. the ODE with respect to y. We then proceed to solving it using a matched 
asymptotic procedure first described in [S] , with some modifications, in section 14.21 In 
addition, we extend the analysis to next-to-leading order. The result of the analysis is a 
relation between the separation parameter A and the parameters m^, and along 
with a new integer n. In section 14.31 we apply the Continued Fraction Method to the 
radial equation, resulting in a discrete spectrum for the Quasinormal modes, which we 
compute numerically. 



The Angular Equation 
The equation, depending only on the angular coordinate y, is 

\d, (vYd,S(y)) - Uv'(-'-*~-)^-'-* + 'M' + 2^ + A I S(y) = 0. (46) 



^2 „,2 



If one uses the variable u = ^-^r^, the equation has four regular singular points. Thus, 
it can be brought to the form of a Heun equation 

r{u) + (I + ^ + ^) F'iu) + F(.) = 0. (47) 

\u u—1 u — sj U[U — S)[U — 1) 

with the four regular singular points located at r = 0, 1, s, oo, and e = {an + Ph — 1 ^ 
5 + 1). The transformation to the canonical form of Heun's equation is performed using 
the following ansatz 

S{u) = {u- lyvP (^u + ^ - 1^ F{u) (48) 

where the newly added constants are chosen so that they cancel the higher order poles 
in the coefficient of F{u). Without loss of generality, we can use the following values 
for them: 

. = ^ . = ^ (49) 
The additional regular singular point, s, is given in terms of a by 

s = 1 - 4 < 0. (50) 
The other parameters of the Heun equation are given by 

a// = 2 + /i + ?7 + X = ^ (4+ |m^| + |m^| - w) (51) 

= lJ' + 'n + x = ]j^i\'m4>\ + \m^\-uj) (52) 
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7= l + 2r/ = l + |m^| (53) 

6= l + 2^ = l + |m^| (54) 

e = 1 + 2x = 1 - w (55) 

q = ^^^^ + ^ (-2u; + 4 |m^| + 2 - uj\m^\ + 2\m^\ + 2\m^\\m^\ + m|) + 



(A + 2 \m^\ + 2\m^\ + 2\m^m^\ + m| + m^) 



(56) 



4a2 

The physical region y G [0,a] corresponds to the canonical interval u G [0,1], with 
y = corresponding to -u = 1 and y = a to u = 0. The point y = I corresponds to 
u = s < 0. In the limit s — )■ two singular regular points merge, and the equation 
becomes hypergeometric. As a result, we can use an asymptotic matching technique to 
solve the equation. 

4.2. Matched Asymptotic Expansion for the Angular Equation 

We will follow the method of [H], refined in [2]. However, some modifications are required, 
since the Heun parameters here are no longer constant, but rather depend here on an 
unknown quantity u. They also contain the parameter |s| itself, and thus must also 
be included in the perturbative expansion. The asymptotic matching is performed by 
dividing the interval u G [0, 1] into two overlapping regions, the far region, where m 3> \s\, 
and the near region, with m <^ 1. The equation in each of these regions reduces to a 
Hypergeometric equation, and is thus solvable. The overlap region u ~ is then 

used to match the two solutions, obtaining a quantization condition. A similar method 
was also used in [33l |3ll |35l [361 EZ] to obtain analytic expressions for quasinormal modes 
in flat backgrounds, where the relevant equation is a Confluent Heun Equation. 

The far region -u 3> |s| Here we can set s = in fHTj) to obtain the zeroth order 
equation, this amounts to setting a = 1 which, in terms of the Heun parameters, will 
affect only q: 

qo = q{a = 1) = r/(l + ri) - + x{l + 2ri + x) - \x - ^{uj - m^) \ (57) 
After some algebra, the solution in this region is given by 

Ffariu) = uf"Fi[aH + Pl, + Pu S, 1 - u]. (58) 
Here Fi is the Hypergeometric function and 



^' = 2 
_ 1 

~ 2 
1 

P2 = - 



an 


-Ph 


+ 6 + 


A/-4go 


+ ( 


-an - 


Ph 




1 - 




-2x- 


f y/l + 


ml 


+ A + ( 


[u- 




an 




+ 5- 


\/-4go 


+ ( 


-an - 




+ 6y) 



(59) 



(60) 
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satisfy the relation 



Pi + P2 = S - an - I^H- 



(61) 



We have also used the freedom to normalize the constant coefficient of this solution to 
1 relative to the near region solution. 

The near region m <^ 1 In the near region, in order to take the limit s — )■ correctly, 
we need to first make the transformation u = in ( l47l) . Then, taking the limit s — > 0, 
one obtains the zeroth order equation in this region. The solution to this equation is 
given by 



where C is a normalization constant to be fixed by the matching procedure. 

Matching the two solutions The next step is to match the two solutions at the point 
u = yj—st with if: ~ 0(1), expanding around s = 0. This is done using hypergeometric 
functional identities. We have already performed the general procedure in [2J, and 
obtained the following quantization condition 



In order to extract the correct quantization condition for A from this equation we must check 
several cases of which only one leads to a true solution. Note that the elimination of other 
possibilities is not completely generic as there is a point in our derivation that depends on 
the values of the parameters of our specific equation. We describe how to eliminate the other 
possibilities in Appendix A, here we focus only on the relevant case. 

Assuming Re(/9i — p2) > 0, the RHS of (j63|) goes to zero as s — )• so the denominator of the 
LHS must diverge in this limit. This means that one of the arguments of the gamma functions 
in the denominator should be equal to a negative integer. These are given, in our equation, 




(62) 



[pi - p2] r [7 + p2] r [-pi] r [/3h + P2] r [aH + P2] 
r [an + pi] r [Ph + pi] r2 [p2 - pi] r [7 + pi] r [-ps] 



2pi-2p2 



(63) 



by 




(64) 



(65) 



(66) 



(67) 



(68) 
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Checking these terms one by one we observe that only the three last terms can be equal to 
a negative integer. In particular, focusing on the third and forth term and assuming that 
Re a; > for simplicity, we get that the eigenvalues are 

A = 4n(l + n) - (2 + 4n) (w - \m^\) - m^^, (69) 

subject to the restriction 

Rew > 1 + Im^l +2n. (70) 

Note that since the condition for superradiance is Re a; < m^Jl, the restriction ()70p means 
that indeed, superradiance is possible only for $7 > 1 and only for m^s satisfying 

2n + 1 

> (71) 



There is a second more technical restriction, in addition to (j70p . This comes from the fact 
that (j63p also gives a higher order (non analytic) correction to the zeroth order eigenvalue we 
found. In particular, taking 7 + = — n + Sn and plugging this into the equation, we get the 
following: 6n oc s''^~p^ which implies that the correction to the eigenvalue itself will also be of 
this order. In order to insure the validity of a perturbative expansion in s, up to first order, 
where this correction may be safely ignored, we require that Re(/02 — Pi) ^ 2. This means that 
eigenvalues with 

Rew > 3 + Im^l + 2n (72) 

are completely valid. Taking into account the non analytical correction will make the entire 
computation harder since the expression p2 — pi depends on u. In this paper we perform the 
full computation without the correction, and discard those eigenvalues that do not respect this 
bound. 

Subleading Correction We would like to go beyond the zeroth order in s = 1 — ^ in 
calculating the angular eigenvalue since in the zeroth order il^ = 1, and we thus expect 
superradiance to be apparent only in the next-to-leading order. To calculate the first order 
correction to the eigenvalue one needs to take into consideration the next order of the 
differential equations appearing in the two regions. The correction can be found separately 
in each region using perturbation theory and a comparison between them is needed to ensure 
this is a valid solution in the entire segment. We describe this calculation in both regions in 
[Appendix B and we only cite here the result (assuming > ): 

Ai = -2 + 2n + 2n^ + (1 + 2n){m^ - uj) + 

{{uj - m^) + (1 + 2n) {u - va^) - 2n(l + n)) (m^^ - 4) 
(2n - cj + m^) {2 + 2n-uj + m^) 

Notice that the last two terms, which are the only ones non- linear in w, cancel when = 2. 

We have implicitly assumed in taking the limit s — t- that all other parameters, and 
specifically m^^ are 0(5*^). Since we have = 1 + 0(s) when a > 2, superradiance can occur 
only for ~ 0(l/s), as superradiant modes must also satisfy (fTTI) . Thus, superradiance, 
and the instability associated with it, appear to lie outside of the domain of validity of this 
approximation. We discuss this issue in greater detail in section 15.11 
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4.3. The Radial Equation 

The radial equation, depending only on the coordinate x, is given by 

1« ^ T-^j ^^ , ((a;2 + 2msf) (-L^ + am^) + am^-a2a;)^ a^m^ 

-dx {xXd^R{x)) + ^ f- 

X \ X x'^ 



Recall, from section [321 t^i^-t upon making the coordinate transformation z = x^, the inner 
event horizon is located at z = while the outer event horizon is at Z2- These, along with the 
boundary of AdS at z — )• 00, are three of the equation's regular singular points, the last one 
being located at 



4a + (a^ + 1 + 22) a 
zi = ^ '—. 75 

a 

Using the transformation r = we move these to the canonical set of points {0, 1, s, 00}, 

with 

. = 1-^. (76) 

^1 

One can now use the ansatz 

R[r) = r^r - (r - I + F{r) (77) 



with 

u = 2 (78) 
_ i [-0J (2a + a (a^ + Z2)) + am^ (2a + q (1 + Z2))) 

2a^{-zi + Z2) ^ ' 

P=:^P^ (80) 

to transform the equation to the canonical form of Heun's equation (j47|) . The Heun parameters 
are then given by 



an = fJ. + V + p - i\l -—fj. (81) 

Zi 



Ph = P + i^ + P + i^-^^f^ (82) 

7 = - 1 + 2i/ = 3 (83) 

6= 1 + 2/i (84) 

e = l + 2p (85) 

(6ia (2 + aa))uj + 6^aam^^/^^Zl + a'^am'i^/z2 + 



Aazi^/z^ 

a (a + + 16zi + 6iuJ^/z^ - 8^2) Vz2- 

2iam^ (3/(2a + a) — iaoj^fz^ + 3az2)) • (86) 

In order to solve this equation numerically, we use the Continued Fraction Method [38], 
which has first been applied to black hole physics in [391 HO] • We first plug into the equation 
a series solution, F(r) = X^^o'^'^'^"' satisfying the regularity condition at the boundary of 
AdS, located at r = 0. This yields a three-term recurrence relation for a„ 

/fcOfc+i + gkO-k + hkak+i = 0. (87) 
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The coefficients are 

fk = sik + l){k + j) (88) 

9k = -k[s{k + -/ + 6-l) + k + -/ + e-l]-q (89) 

hk = {k + aH-l)ik + f3H-l). (90) 

We also require regular boundary condition at the other end of the physical interval, r = 1, 
which corresponds to the event horizon. This requirement is only satisfied by the minimal 
solution which in turn requires convergence of the following continued fraction 



n /o • /ll /o • hi fi • /l2 /2 • /l3 ^ 

= 90 — = 90 • (91) 

/i-/i2 91- 92- 93- 



91 



f2 ■ hs 

92 



93 

This sets the quantization condition for u. 

Truncating the continued fraction after some p steps, we get a polynomial equation for 
uj which can be solved numerically. Increasing the cutoff causes more quasi normal modes to 
be found, as well as better accuracy for the ones already found. This is true as long as the 
precision of the computations is increased as well, so that numerical errors are minimized. 



5. Scalar Quasinormal Modes 

Using Mathematica we carried out the numerical procedure for various choices of the 
parameters. In order to achieve numerical stability, a large cutoff of the continued fraction 
was used, often over 140 terms, in high precision floating point arithmetics. Below we present 
an example of the typical results obtained using the numerical procedure. 

As expected, for small m^s, we did not find any instabilities using a variety of values for 
a, m^, m^, Z2 and re. As an example. Figure [1] displays the first quasinormal modes for the 
following parameters a = 2, a = 0.99, Z2 = 0.2, = 7, = 2, n = 1. 



3 -15- 



10 15 20 25 30 35 

Re[w] 



Figure 1. First 8 QNMs for a = 2, a = 0.99, Z2 = 0.2, n = 1, = 7, = 2. The 
continuous blue line denotes the energy bound ([70]) . the dashed red line denotes the 
bound for perturbation validity (|72p 
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The relative numerical error resulting from the cutoff of the continued fraction is less 
than 0.0001. 

5.1. Observing Superradiant Instabilities 

In order to observe superradiant modes, we require uj < m^fl^, but we also know that the 
approximate solution is valid only for to > + 1. Thus, we must have rrifj, > oc jiy. 

Having such a large parameter in the problem turns the angular equation into the more 
complicated Confluent Heun Equation. However, naive computations using the numerical 
method indeed show the known superradiant instabilities, as found in other types of black 
holes O El [HI [171 HH US] As an example of this superradiant instability, we display the first 
quasinormal modes for the following parameters a = 2.1, a = 0.9999, Z2 = 0.01 ,n = 0, 
= lO'^, = 2. In this case, we have fi^ — 1 ~ 1.9 x 10^^. We can estimate the errors to 
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Figure 2. QNMs for a = 2.1, a = 0.9999, Z2 = 0.01 ,n = 0, = lO"^, = 2. The 
continuous blue line denotes the energy bound (|70)) and the dashed red line denotes 
the bound for perturbation validity (|72|) . The thick black line shows the expected 
superradiance threshold. 

this numerical computation by increasing the cutoff on the continued fraction, and checking 
the change in the eigenfrequencies. In this case, varying the cutoff from p = 100 to p = 140 
shows that the real part of the frequencies is correct to about 10~^, while their imaginary part 
is less accurate. Far enough from the crossing point between positive and negative imaginary 
part we have an error of approximately 0.05, while close enough to the superradiance threshold 
(depicted by the black line in Figure [2|), errors grow up to about 1. This explains why the 
actual crossing point is not precisely at the superradiance threshold. 

The well-known phenomenon of superradiance thus occurs where expected. The reason 
that the scaling ~ does not mask it, is that even with this scaling, we still expect 
(to — m^) to be finite without introducing new singularities. Indeed, if we search for eigenmodes 
having a; = + /S. with A ~ 0(1), Ai is still finite, with the following exception. There is one 
term that, for 7^ 2, is of order 0(l/|s|) (the term is absent for = 2, which is depicted in 
Figure[2]). This term changes the zeroth order term in the expansion of A, but it remains finite. 
Thus, the result for A, obtained from the next-to-leading order of the asymptotic matching 
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procedure, is finite in the limit s — )• 0. If tliere are no further corrections to the zeroth order 
term coming from higher order terms when using this scaling, the results obtained here would 
be correct at leading order. We have not proven that this happens, but the compatability of 
our results with the expected superradiance threshold support it. 



6. Thermodynamic Stability 



Hmn = - ; = {E,J^ = J,Q) (92) 



a = — (93) 



As previously mentioned, thermodynamic instabilities manifest themselves as non-positive 
modes of the thermodynamic hessian 

dx'^dx^ 

In order to compute the Hessian, we first express the entropy using these variables, by inverting 
(jl6p . Then, after computing the Hessian, we would like to express the result using the 
parameters a, a and Z2. However, it is more insightful to express the Hessian in terms of 
a,6,Q where a and 5 are defined in the following way 
J_ 

and make a large-charge expansion, taking the limit Q — )■ oo. This limit is equivalent to the 
limit a — >• 1 we have considered thus far and in taking it we must keep a,5 ~ 0((5'^). At 
leading order in Q it is easily verified that the extremality is achieved for b taking the value 
Z = We thus substitute 6 = S + Sq, so that states above extremality correspond to values 
of (5o > 0. The Hessian at leading order is given by 

(a —a —2y/2a \ 

-a a 2V^ • (95) 

The eigenvalues of the hessian are 

Ai = (96) 



4 + a + 2V2a + ^J{'i + a)'^- 4\/2^(a - 4)5o + 8a6^ 
As = ■ ; (97) 



4 + a + 2v^ - y (4 + a)2 - 4^/2^(a - 4)(5o + 8a5'^ 

As = ■ ; Tj:^ (98) 

25/4a9/4g53/s 

Thus, at leading order the hessian always has a zero mode, and we must go to next-to-leading 
order in Q. Since we know As and A3 to leading order, the determinant 

allows us to extract the leading order contribution to Ai, and we find: 
_ 2-a + V2^5o _ 

2l5/4Q3«9/4^ 2lg2„3/4^ ^ 
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4 + a + 2V2a - a/ (4 + a)2 - 4^/2^(a - i)5o + 8a6?, 

M = TT. > (102) 

25/4a9/4g53/2 

where when expressmg Ai in terms of we used the expansion of Q.^ up to order O(^). 
The result (jlOOp reflects the fact that the black hole is thermodynamically unstable when 
> 1. This is the same regime where superradiant instability has been shown to occur 
in the previous section, a fact which might be related to a Gubser-Mitra-like conjecture [TU] 
for rotating black holes, though there are some differences, as discussed in the next section. 
The eigenvector corresponding to this instability is given by (1,1,0), when written in terms of 
{5E,6J,6Q). 

In order to describe the onset of the instability in terms of the charges E,J,Q, let us use 
them to re-express Q^: 

-E + V2y/J + J / 1 

Thus, the stability would appear for E — J < 2\/j. In particular, for very large energies there 
is no instability while for very large angular momentum instability is guaranteed. 
It is also interesting to examine the temperature in terms of the charges 



^0 = 1 + +(7^)- (103) 



J {-2E^/J + V2 J + 2 J3/2 + 2V2Q2 



2i/4j^ 

Therefore, as we could have also concluded from (j95p . the temperature is an increasing function 
of energy at fixed J and a decreasing function of the angular momentum at fixed E. 
For completeness, we express 5 in terms of the parameters a, a, Z2 

4a2(3 + a) + 4a (2 + a + a2 + (3 + a)z2) a + (3 + a) (1 + Z2) (a^ + Z2) o? 



2^/2(1 + a)^Jaa (2a + q + Z2a) {2a + (o^ + Z2) a) 

I +^^ + 0(1 -a) (105) 
v2 



7. Summary and Discussion 

We have analyzed the stability of a three parameter family of black holes charged under two 
of the U{l)s in the cartan of 50(6), within the consistent truncation of type IIB SUGRA on 
AdS^ X S^. The specific black holes that we studied had two equal charges (denoted by Q). 
The black holes also have non-zero angular momentum around one of the two independent 
rotation planes of AdS^, denoted by J^. We have defined a = N'^J^/Q'^ - a quantity measuring 
the ratio between the angular momentum and the charge squared. For this class of black holes, 
as we approach the extremal limit, the horizon shrinks to zero size. Such black holes have 
been dubbed Extremal Vanishing Horizon (EVH), and their properties have been studied in 
[HEIIHD]. For a generic point on the horizon, the extremal near horizon limit is pinched AdS^. 
In contrast, on the plane of rotation the horizon merges with a singularity as we take the 
extremal limit, and no universal near horizon exists at the level of the geometry. The value 
a = 2 is distinguished in the sense that it satisfies, at extremality, a SUSY BPS bound. 

We carried out a stability analysis in the limit of J,^ — )• 00 for fixed a. In the notation of 
the solution given in this limit corresponds to the limit in which the rotation parameter 
a — )• 1. We have checked for both local thermodynamic instabilities and classical instability, 
analyzed using linear perturbations analysis. The classical stability analysis was carried out 
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for a scalar field, after identifying an appropriate scalar whose addition still forms a consistent 
truncation of type II B SUGRA on AdS^ x S^. 

We found evidence for both types of instabilities when > 1, where fi^ is the 
angular velocity of the black hole with respect to a non-rotating frame at the boundary of 
AdS. This is also the threshold predicted in [9] for superradiant instability, induced by the 
asymptotically AdS nature of spacetime. Indeed, the unstable modes are precisely those which 
are superradiant. In the near extremal limit this condition translates to a > 2. 

It is important to note that, in the limit a — t- 1, the angular velocity takes the form 
0,^ = 1 + 0(1 — a) and one must therefore examine the sub-leading term to see if > 1 or 
not. In addition this form of 17,^ implies that a superradiant instability can appear only for 
modes of the scalar having very large angular momentum, m,^ ~ 1/(1 — a). 

It would be interesting to find out what would be the endpoint of these instabilities. 
Unfortunately, knowing there is a thermodynamic instability does not allow us to extract much 
information regarding the final state, and the classical linear stability analysis is obviously 
insufficient for such purposes, and should be replaced by a full non-linear analysis. 

A natural question to ask at this point is whether the two instabilities are related. It has 
been shown numerically that for AdS^ Reissner Nordstrom (RN) black holes, in the brane limit, 
local thermodynamic stability coincides with dynamical stability |43t I1U|. This phenomenon 
was shown to persist also for black holes with horizon radius comparable to AdS radius. The 
authors further suggested that it is a generic phenomenon, in what came to be known as 
the Gubser-Mitra conjecture [10]. A part of this conjecture, which states that black branes 
are unstable if, and only if, they are constructed from a thermodynamically unstable black 
hole (at lower dimension) , has recently been proven in |34J . This work also established a close 
connection between dynamical and thermodynamical stability of black holes, but unfortunately 
it has not yet been extended to Anti-de Sitter space. It also does not apply to our case, since 
the superradiance instabilities that we find are at short wave length in the J — t- oo limit. 

In the context of fast-rotating black holes, which we are discussing here, intuition can 
be obtained from arguments originating in the study of higher dimensional asymptotically 
flat black holes. It has been argued [45j that black holes in L> > 6 with one non-vanishing 
rotation parameter are prone to ultra-spinning instabilities. Black holes are ultra-spinning 
when there is no upper bound on the angular momentum for fixed mass. For large enough 
angular momentum, the event horizon becomes quasi-extended, resembling a black brane. 
The observed instabilities, which are purely gravitational, are in fact Gregory-Laflamme (GL) 
type instabilities causing fragmentation (or clumping) of the extended black hole. As 
mentioned above, for black branes, such instabilities can be identified by studying the local 
thermodynamic instabilities. This can be understood intuitively by arguing that the clumped 
black holes are entropically favored when a local thermodynamic instability exists. 

As mentioned in [35], the above argument cannot be made precise m D = 5 since the 
black holes in that case have no ultra-spinning regime, and the black brane limit is not 
defined since there are no black two-branes. However, fast-rotating black holes in D = 5 
asymptotically flat spacetime can still become very much 'pancaked', as is the case for the 
black holes studied in this paper. It has been conjectured j4_7J that such black holes become 
unstable before extremality, with the endpoint of the instability being a black ring. This 
black ring is thermodynamically favored, so that a thermodynamic instability is expected 
for the black hole. In [38] a classical gravitational non-axisymmetric instability was found 
which approximately matches the instability region predicted by thermodynamic arguments. 
Whether this instability also occurs in AdS^ black holes is, to the best of our knowledge, 
unknown. We have not yet performed a gravitational stability analysis, so it is not clear if 
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such an instability exists, but it is obviously a viable option. 

However, this is not the only possible outcome for the black hole evolution. Superradiant 
instability, which is mostly absent from the discussions in the above cases, occurs also at zero 
temperature (i.e. in extremal black holes). It can be interpreted within string microscopic 
models of black holes as a quantum effect [49J. At zero temperature, superradiance is taken to 
be spontaneous emission, while at non-zero temperature, the thermal radiation is thrown into 
the mix. It is possible that the instability observed in the thermodynamic analysis reflects this 
as well. The stable configuration, if this is the dominant mechanism of decay, is most likely of 
the type predicted in PT] for four dimensional AdS black holes. It is a slightly modified il^ = 1 
black hole co-existing with some outside scalar hair and radiation. A scalar field instability 
of small charged and rotating black holes in AdS^ was shown numerically to result in a black 
hole with scalar hair ^50j (also [51j). Hairy black holes were also shown to occur in AdS^ x 
[52] . In our case, we expect that the final result would be an a = 2 core surrounded by hair. 
Of course the analysis for the black holes under consideration here is further complicated by 
the fact that they have a non-trivial configuration of both the gauge fields and scalars. 

The black holes which are the focus of this paper have been conjectured in p] to 
correspond to Fermi surface-like operators in the PS'C/(1,1|2) sector of the dual conformal 
field theory. A/" = 4 SU{N) SYM. Understanding the final state evolving from the unstable 
black holes can help one study the dual CFT vacuum configuration for a > 2. Conversely, if 
this instability can be observed on the field theory side, perhaps it would hint at the end-point 
of the gravity theory's instability. 

The classical instabilities observed in this paper are for neutral fields, carrying no charges 
out of the black holes. It would be interesting to see if the black hole can also emit some of its 
charge to the outside. This can be studied by either adding charged scalars to the consistent 
truncation, or using charged vector fields, using the truncation to SU{2) x U{1) instead of 
C/(l)^ [S3]- Another possibility would be to compute Fermionic two-point functions, which, in 
addition to charge emission, can also be used to observe Fermi surfaces in the dual CFT. Such 
fermi surfaces have been observed for the black-brane analogs of these black holes [54] . 
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Appendix A. Elimination of additional possible quantization conditions for 

A 

We need to consider whether the relation p2 — pi = —n + 0{s) might lead to some quantization 
condition. This was treated rather generally in the Appendix of our previous work, up to a 
point which depended on the parameters of the equation. We will treat this point here, using 
the notations of that Appendix. The non generic case to consider is when a = —m or b = —m 
and n < m for both the far and the near region in which case the expansions in s leading to 
()63p are valid, a and b for the far region are defined by 



afar = f^H + P2 

bfar = OH + P2 (A.l) 

while for the near region they are 

Oinear — Pi 

bnear= - Pi - J + I = -\in^\ - pi (A. 2) 

Now we can rewrite (l63l) as 



[n] r [7 + P2] r [amar] r [afar] T [fe/gr] ^ .2n . . oN 

T[aH + Pi]T[f3H + Pi]T^-n]T[j + pi]T[-p2] ^ ^ 
We observe that for this equality to be fulfilled we need to have at least one more divergent 
gamma function in the denominator than in the numerator of the LHS since the RHS goes 
to zero. We will show this doesn't happen, assuming throughout our considerations that the 
only divergent term in the denominator is (F [— n]) ^ since the other cases were already treated. 
Thus, we need to check only the cases where there is only one divergent gamma function in 
the numerator 

The only such case is that one of afar, bfar is equal to —mi,—m2 > —n and bnear also 
fulfills this requirement but a„ear doesn't fulfill it. This means (1) — — pi < and (2) 

+ Pi| = w-i < ^- From (2) pi is an integer. If —pi < then from (2) < n so flnear 
fulfills the conditions contrary to our assumptions. If pi = then F[— pi] is divergent so the 
LHS is finite and incompatible to the RHS. Therefore pi < and from (1) —m^ ~^\pi\ < =^ 

< m<^. Now we can rewrite (2) as — |pi| < n so m,^ < n. Recalling that pi = n + p2 
we get that p2 is an integer and using the fact that pi < we get that p2 < —n. Now we shall 
prove that the in this case F[7 + P2] is divergent. 

7 + P2 = H-m0 + p2<H-m^-n<l (A. 4) 

since both 7 and p2 are integers we got that 7 + P2 is an integer at most equal to zero and 
this ends our proof. 

Appendix B. First order correction to the A in the a — )■ / case 

The far region 

In the zeroth order we had the equation 

Dotpniu) = qo{n)^n{u) (B.l) 
with the differential operator 

D^ = u [{u - l)udl + {{u - 1)(1 + + (5h) + 5)8^ + (B.2) 
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and the solution was given by 

V'" = uf"Fi[aH + pi,pH + Pi,S,l-y] 

= u'^-^'Fi [an - - n, pH - 1 - n,5,l - u] (B.3 

after using the quantization condition from the matching. 

Thus, taking F{u) = '0° + as the solution one gets the equation up to first order in , 
(Do + + si^l) = (qoin) + sq,){i^l + (B.4 

which gives the following equation for the unknowns -0^ and gi 

D^i^l - q^{n)i^l = qi€ - (B.5 
One can actually rewrite the differential operator of the full equation as Dq + cDi where 

= u{l - u)dl + ((1 + aH + Ph){1 -u)- 6)du - e(l - r)d^ (B.6 

^--DQ-e{l-u)Dr + aHPH (B.7 
u 



Then, using the following recurrence relations 



with 



kn+1 = 7 w TT (B.9 

Pi - P2) [Pi - P2 - 1) 



-I + /3h + Pi + P2 - 2piP2 - ;0JJ (Pl + P2) 

(Pl - P2 - 1) (Pl - P2 + 1) 
Qijj (-1 + 2/3g + pi + P2) 

(Pl - P2 - 1) (Pl - P2 + 1) 

(aH + Pi)(/3g + Pi) 

(Pl - P2) (Pl - P2 + 1) 



and 



i=-l 



with 



(B.IO 



K-i = + + (B.ll 



(U - l)dMu) = E Pn+^V'°+^ (B.12 



(ag + Pi)(/3g + Pi)p2 .p. TO 

Pn-l = -7 T-TT- T (B-l<3, 

(Pl -P2) (1 +Pl -P2) 

P1P2 (Pi + P2 + -l) + aH (2pip2 + (1 + Pl + P2)) 1 , 
(pi - P2 - 1) (pi - P2 + 1) 
{aH + P2) + P2) Pl T^ 
(pi - P2) (-1 + Pl - P2) 

So finally we can write 

i=+i 

{Do - qo{n))i^l = Qii^l - ^ i^Pn+i - kn+iqo{n))i^l^i (B.16) 

i=-l 
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Wc must use a solution of the form i/jI^ = An+iU^,^.^. Since -0° are eigenfunctions 

of Dq with different eigcnvahics they are hnearly independent and we get an algebraic 
equation for the coefficients Aj. Examining the coefficients of m„ we get an expression 
for qi 

qi = aH/^H + epn - qQ{n)kn (B.17) 

This means we got the next order of the eigenvalue without using the matching of the 
two regions. To make sure this makes sense we should perform the same procedure in 
the near region and compare the two corrections. 

The near region 

The zeroth order equation is 

DoXniO = Qoin)XniO 
with the differential operator 

Do = + Od^' -{l+{l + aH + /3H- S)Od^ 
and the solution was given by 

= Fi (7 + n, e - 1 - n, 7; ■ 
The next order equation is 

Doxli - (lo{n)xl = qixl - Dixl 

with 

Di = + 0^?' - e(7 + <^ + (1 + «H + MOD^ - o^hPhC 



One can express 

CXn ^ ] fn+iX. 



i=+l 



i=-l 



with 



fn-l 
fn 
fn+1 



{pl + 7) P2 



ip2 - Pl - 1) {P2 - Pl) 
2piP2 + 7 (1 + Pl + P2) 
{p2 - Pl + 1) (P2 - Pl - 1) 

- (P2 + ^/ ) Pl 



and 



B.18) 
B.19) 

B.20) 
B.21) 

B.22) 
B.23) 

B.24) 
B.25) 
B.26) 



(P2 - Pl + 1) (P2 - Pl) 

i=+i 

+ l)Dan = ^n+iXn+i (B.27) 



i=-l 
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with 



=n-i = -P^P^i^ + P^) ^ (B.28) 



in+l 



{P2 - Pi) (P2 - Pi - 1) 
P1P2 (P1+P2 + 27-I) 
(P2 - Pi - 1) (P2 - Pi + 1) 
-P1P2 (7 + P2) 



(B.29) 
(B.30) 



(P2 - Pi) (P2 - Pi + 1) 

Using these expressions and writing Xn as a hnear combination of Xn+i ' u^ith — l<i<l, 
in (]B.21I) . one can proceed in the same manner as in the far region and obtain 

qi = (qoin) - aHPH)fn - 5e„. (B.31) 

Indeed, comparing the two expressions for qi from the far and near region, using 
qo = P1P2, one gets an equahty. We can now proceed to compute the first order correction 
for A from the expression for qi . We have q = qo + sqi with 

go = ^ (2^0 - - Ao + 2w + im^u) (B.32) 

and 

gi = ^ (2m0 + m^^ + 2m^ + 2m^m^ + m^^ + Aq - m^u) - ^Ai. (B.33) 
Plugging this into f lB.3ip and solving for Ai we are left with 

Ai = -2 + 2n + 2n^ - (1 + 2n)u + + 2nm^ + 

(-4 + m^^) n (n + m^) (-4 + m^^) (1 + n) (1 + n + m^) 



(B.34) 



2n — u + niff, 2 + 2n — u + 
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